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We show that in a macroscopic perfect insulator, charge 
injection at a field-enhancing defect is associated with an in- 
stability of the insulating state or with bistability of the insu- 
lating and the charged state. The effect of a nonlinear carrier 
mobility is emphasized. The formation of the charged state 
is governed by two different processes with clearly separated 
time scales. First, due to a fast growth of a charge-injection 
mode, a localized charge cloud forms near the injecting de- 
fect (or contact). Charge injection stops when the field en- 
hancement is screened below criticality. Secondly, the charge 
slowly redistributes in the bulk. The linear instability mech- 
anism and the final charged steady state are discussed for 
a simple model and for cylindrical and spherical geometries. 
The theory explains an experimentally observed increase of 
the critical electric field with decreasing size of the injecting 
contact. Numerical results are presented for dc and ac biased 
insulators. 

PACS number: 77.22. Jp 



I. INTRODUCTION 

Insulation of dielectrics is limited due to dielectric 
breakdown There exist several different physical 

mechanisms which lead to instabilities associated with 
dielectric breakdown at high electric fields, e.g., thermal 
runaway and impact ionization avalanches. At break- 
down, a change from an insulating to a conducting state 
occurs, at least in a certain spatial region and for a cer- 
tain time. A release of charge carriers is possible from 
two different sources. First, carriers can be generated 
intrinsically via a bulk instability, e.g., by ionization of 
impurities. Secondly, carriers can enter due to injection 
at the electrodes. In this paper, we show that charge 
injection is associated with an instability, too. In con- 
trast to bulk instabilities, however, charge injection is a 
boundary instability where the unstable mode (charge 
injection mode) is localized at the injecting contact. 
In practice, charge injection in macroscopic insulating 
bodies occurs at geometrical defects of the electrodes 
where the electric field can be strongly enhanced. Be- 
low we will consider concentric cylindrical and spherical 
contact geometries. A small inner electrode has a large 
electric field and can serve as a model for a field enhanc- 
ing defect. The cylindrical system describes also a coaxial 
cable filled with a dielectric medium; a system which is 
of obvious technical interest. 

Charge injection in dielectrics has been investigated ex- 



perimentally for a tip-plate geometry by Hibma and 
Zeller Q], and has been modeled by Zeller and Schnei- 
der g| in the limit of an infinitely high mobility edge 
and by neglecting diffusion. They describe the mobility 
edge with a mobility [i(E) which vanishes for E < E c 
but which is very large for E > E c , where E c is a critical 
value of the electric field (~ 10 7 V/cm). In this model 
a space charge forms near the tip when the local field 
reaches the mobility edge. The space charge, in turn, 
screens the electric field enhancement at the tip and pins 
it to the mobility edge (field limiting space charge). Zeller 
and collaborators assume a bulk instability at E c asso- 
ciated with S-shaped negative differential conductivity 
which forms the basis of their theory of injection [g||3|. 
Below we show, however, that there is no need for such 
an underlying bulk instability. Charge injection turns 
out to be an instability by itself. 

Boggs |?],|| pointed out the usefulness of the screening 
by the injected space charge in ac driven field-grading 
materials. His model, however, is based on the concept 
of conductivity which cannot lead to a consistent physi- 
cal description of charge injection. The theory assumes 
a conductivity which is only a function of the field and 
which does not distinguish between intrinsic and injected 
charge carriers. Note that charge injection is a boundary 
effect, while conductivity is a bulk quantity associated 
with intrinsic carriers. Boggs' approach leads, nonethe- 
less, to qualitatively correct ac results in the limit of an 
infinitely sharp mobility edge and in a certain frequency 
regime. 

For the sake of clearness, we will consider a perfect insula- 
tor, which is defined here as a dielectric without intrinsic 
carriers and with a constant permittivity e. Electrons or 
holes can be present only due to injection at the elec- 
trodes. Without charge injection, a voltage difference at 
the contacts induces a charge which is located outside 
the dielectric in a thin surface layer (with a thickness of 
the very short Debye length) of the metal contacts. We 
call this locally neutral state of the insulator the ideal 
insulating state. The electric field in the insulator is then 
uniquely determined by the Laplace equation. The elec- 
tric field at the contacts is fully determined by the po- 
tential differences of the contacts. Clearly, this situation 
corresponds to a purely capacitive arrangement of the 
contacts in a dielectric medium. 

A prescription of arbitrary boundary conditions to the 
electric field at the contacts (e.g., E = 0) is more restric- 
tive and in general implies the presence of a space charge 
in the dielectric medium. The field is then determined 
by the Poisson equation. Important work on charge in- 
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jection ]l|,p|-|T2]| treats the formation of a space charge in 
this way as a direct consequence of boundary conditions. 
These theories treat the charged state, but they do not 
consider the stability of the locally neutral state. 
A different approach which is appropriate for metal- 
semiconductor contacts is to prescribe a Richardson- 
Schottky @y or a Fowler-Nordheim Jl5|]l6| ] current- 
field characteristic in order to model thermionic (field) 
emission or a tunneling current through the contact bar- 
rier, respectively. In contrast to the well-defined metal- 
semiconductor micro-contacts manufactured by a highly 
developed semiconductor technology, macroscopic metal- 
insulator contacts used in high-voltage devices are not 
well-defined and can thus not be treated on a microscopic 
level. A description on a hydrodynamic level is then more 
appropriate. Therefore, we prescribe phenomenological 
boundary conditions to the charge density p. Additional 
boundary conditions to the electric field are unnecessary. 
In principle, the parameters occurring in the boundary 
conditions should be derived from microscopic models, 
provided the physics at the contacts is known. For ho- 
mogeneous boundary conditions, it turns out that the 
locally neutral state (p = 0) is always a stationary solu- 
tion of the problem. However, we will show that this ideal 
insulating state can become unstable against a charge in- 
jecting mode or that bistability of neutral and charged 
state can occur. 

This paper is organized as follows. In the next section, 
we introduce a model for the perfect insulator with phe- 
nomenological boundary conditions. In Sect. Ill, we in- 
vestigate the charge injection instability of the ideal insu- 
lating state. The steady state which eventually developes 
is discussed in Sect. [V. Finally, in Sect [v|, we present 
numerical results for an ac biased perfect insulator. 



II. THE MODEL 

We consider a material with electron-hole symmetry 
and with an immediate recombination of electrons (n) 
and holes (p) by annihilation. This means that the 
mobilities and the diffusion constants of electrons and 
holes have equal absolute values, p = —p n = fJ-p, and 
D = D n = Dp, respectively. An Einstein relation is not 
considered for the present nonequilibrium system, and we 
assume that D is a field independent constant. Due to 
the fast electron-hole recombination, the carrier density 
is equal to the absolute value of the charge density. The 
drift current is then simply given by p(E)\p\E. Dynamic 
equations for electrons and holes describing generation- 
recombination processes do not appear. We emphasize 
that, except for the ac results, all results below are valid 
also for unipolar conduction and are not consequences 
of the electron-hole symmetry and of the fast recombi- 
nation. In particular, the injection instability discussed 
below occurs for unipolar conduction. 
Consider now a perfect insulator in a capacitor of cylin- 



drical or of spherical symmetry. Metal contacts are at- 
tached at the inner and the outer radius, r\ and ri), 
respectively. In the following, the cylindrical capacitor 
of length L z and the spherical capacitor are labeled with 
d = I and d = 2, respectively. All quantities depend only 
on the radial coordinate, r. The (radial) current density 
can be expressed in terms of the charge density p and the 
(radial) electric field E: 



] = p{E)\p\E - Dd r p 



(1) 



We assume a mobility p{E) which depends on the field 
in the form p(E)E = v\E/ E \ a sign(E), where a > 1 
is a measure of nonlinearity and where v is a positive 
velocity. The limit a — > oo corresponds to the infinitely 
sharp mobility edge at E = E discussed in Ref. || . Note 
that already the case a — 1 corresponds to a nonlinear 
current-field relation since p is related to the electric field 
via the Poisson equation 



eV • E = p 



(2) 



Consequently, a linear dielectric relaxation mode does 
not exist in the perfect insulator. 

There are two equivalent formulations of the dynamics, 
namely in terms of the Maxwell equation 



ed t E = V x H -j 



(3) 



which is a dynamic equation for the electric field, or in 
terms of the continuity equation, 



d t p = -V • j 



(4) 



which is a dynamic equation for the charge density. For 
convenience, we use below Eq. (||) for the numerical sim- 
ulations and Eq. (Q) for the analytical discussion. 
The system is driven electrically via a coupling to an 
external electric circuit, which consists here of a voltage 
bias V{t) and an ohmic resistor, R ex t, in series. The total 
(radial) current density (V x H) r = J/r d is determined 
by 



J = a d [ V(t) 



Edr 



(5) 



where ai = (2TrL z R ext )~ 1 and 0,2 = (47ri? ea; t) _1 for the 
cylindrical and the spherical case, respectively. We men- 
tion that Eq. (||) gives rise to a strong nonlocality which 
can influence qualitatively the spatio-temporal dynam- 
ics of the system p7[ . Below, we restrict ourself to the 
limit case of voltage control, i.e., R ex t — > and to low fre- 
quencies such that inductive effects can be neglected. An 
increase of R ex t corresponds to forcing a current which 
requires the presence of charge and is thus expected to 
lower the stability of the ideal insulating state. Volt- 
age control can equivalcntly be expressed in the form 
V(t) = C iE dr. 

In order to have a well-defined problem we specify mixed 
homogeneous boundary conditions to the charge density 
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d r p\ ri 2 ± Kp\ rii2 







(6) 



where k is a phenomenological parameter, and where + 
and — refers to r\ and r%, respectively. Some remarks 
concerning this boundary condition are in order. First, 
a restriction to homogeneous boundary conditions is not 
necessary. An additional inhomogeneity in Eq. (|^) leads 
to a finite boundary charge. In this paper, however, we 
want to show that charge injection occurs even for ho- 
mogeneous boundary conditions where a locally neutral 
state exists. Secondly, K can depend on the local electric 
field. Such nonlinear boundary conditions can lead to in- 
stabilities. Below we show, that even for the linear case 
an instability occurs, and we discuss the behavior of the 
perfect insulator as a function of k. Thirdly, we assume 
that the charge does not 'wet' the contacts, i.e. k < 0. 
This is reasonable if the microscopic contact potential 
has the shape of a barrier. In a purely diffusive system, 
a 'wetting' density leads to an instability of the uniform 
state. For homogeneous Neumann boundary conditions 
(k = 0) which describe contacts with vanishing diffusion 
current, the ideal insulating state in the diffusive regime 
is marginally stable (gapless stability spectrum) . Indeed, 
an arbitrary spatially uniform p is a solution in the lin- 
ear diffusive regime which implies the existence of a zero 
mode. For finite negative K, the p = state is stable 
in this regime. In the following section we show that, on 
the other hand, an instability of the ideal insulating state 
occurs in the drift dominated regime. 



III. INSTABILITY OF THE IDEAL INSULATING 
STATE 

In this and the following section, we consider a sta- 
tionary and positive bias voltage applied to the contacts, 
V(t) = V > 0. Obviously, a steady state of the system 
is given by p = and E = Cd/r d with V/Ci = ln(ri/r 2 ) 



and V/C 2 



" 1 . This ideal insulating state cor- 



responds to a purely capacitive system. To test the lin- 
ear stability of this state, we seek for the dynamics of a 
weak perturbation (SE,Sp) oc exp(Ai) which satisfies the 
boundary conditions (|B|). From the continuity equation 
(|]) , one finds an eigenvalue equation for the growth rate 
A 

**P + ^ {^fv\S P \ ~ Dr d d r 5p^ = . (7) 

An instability of the ideal insulating state occurs if there 
exists an eigenvalue A with positive real part, since the 
mode Sp associated with such a A grows exponentially in 
time. A dimensional analysis of Eq. (^) leads to a scaling 
relation for the growth rate, 



A = S/(A) 

n 



(8) 



where A = (r\v / D)(E\ / En) a has the meaning of a di- 
mensionless control parameter. Here, E\ — Cd/rf is the 



electric field at the inner contact. Note that the func- 
tion / depends still on d, a, kt\. The dependence on 
KT2 is weak for r\ <C r 2 and will be suppressed whenever 
possible. The critical field at instability depends on the 
various parameters in the form 



E r — En 



D 



-A c (a, d, kt\) 



(9) 



where the function A c has to be determined from / = 0. 
The eigenvalue A with the largest real part turns out to 
be purely real and can be estimated if either the diffusion 
current or the drift current dominates. First, if the drift 
term can be neglected, Eq. (R) reduces to a linear dif- 
fusion equation. Consequently, the eigenfunctions of the 
stability problem are of diffusion type and are damped or 
marginally stable, provided k < 0. On the other hand, if 
the diffusion term can be neglected, the stability problem 
reduces to a first order differential equation. Solving Eq. 
(0) for D = leads to a positive growth rate 
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(10) 



associated with a perturbation 
5p{r) cx r d ^ a -^ exp 



r \ da+1 d(a — 1) + r\K) 



ri 



da + 1 



(11) 



Equation (|llj) describes the unstable injection mode 
which is localized at the inner contact. Obviously, a 
negative k acts to slow down the growth of the unsta- 
ble mode. In Fig. [l] numerical solutions of the stability 
problem of the cylindrical case (d = 1) are shown as a 
function of A with «xi = —0.5. For D/(vri) — > and 
A > 0, the numerical results are in good accordance with 
the approximate analytical results ( |l0| ) and ([ll]) . 
The physical mechanism for the instability can be under- 
stood as a positive feed-back process. Consider a large 
electric field at r = r\ , and assume a negative field fluc- 
tuation, 5E < with d r 5E > —SE/ri, localized at r%. 
The Poisson equation implies then a positive charge fluc- 
tuation 5p at this contact. Using the linearized Maxwell 
equation, edtSE w — Sj < 0, one concludes that the neg- 
ative field fluctuation grows in amplitude, if drift domi- 
nates diffusion. The initial perturbation is thus amplified 
which characterizes an instability. We mention that for 
a = 1, the charge injection mode ([ll]) has no physical 
meaning. In this case, the ideal insulating state is lin- 
early stable, although it is not necessarily globally stable. 
For finite D, a competition between the stabilizing dif- 
fusion term and the destabilizing drift term leads to a 
finite critical value of the control parameter, A c , or equiv- 
alently, to a finite critical field E c . In order to discuss the 
dependence of A c in Eq. (||) on nri, we solve Eq. ([?]) at 
A = 0. A solvability condition leads then to an expres- 
sion for kt\ as a function of A c (appendix). The result is 
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plotted in Fig. ^p, for various values of the nonlinearity 
parameter a, for d — 1, 2, and for constant radii. Clearly, 
A c vanishes for k — + 0. On the other hand, the locally 
neutral state becomes more stable as — k increases due to 
a decrease of the charge density of a density fluctuation 
at the injecting contact. 

In similar way, stability analysis yields the critical field 
E c as a function of r\. We find that the critical field is 
almost independent of r\ except for small n, where E c 
becomes large. This behavior is more pronounced as a 
increases. For a = 3 and d = 1, 2, the results are shown 
in Fig. ||b. For a tip-plate geometry, Hibma and Zeller 
found experimentally that the critical field is almost 
independent of the tip size in a large range but increases 
considerably for very small tip radii. Our theory clearly 
reproduces this behavior. 



IV. THE CHARGED STEADY STATE 

The injection of the charge acts to decrease the field 
enhancement. Consequently, the growth of the injection 
mode saturates at a field below the critical value E c . 
Zeller and Schneider (^] observed that in the infinitely 
sharp mobility-edge limit, a — > oo, the final state con- 
sists of a charged region with p cx 1/r and E(r) w E c 
for r\ < r < f, and a locally neutral region, p = and 
E oc l/r d , for f < r < r 2 . The outer radius f of the 
field limiting space charge is determined by the conti- 
nuity of E(r) at f and by the prescribed voltage drop, 
V = J E dr. One expects for finite a || and in the pres- 
ence of diffusion, that this state decays on a long time 
scale and is in fact part of a transient behavior. More 
concrete, the charged steady state forms on two clearly 
separated time scales. On a fast time scale determined 
by Eq. (||), charge is injected such that the electric field 
drops locally below the critical field. In a second step, the 
charge distributes slowly towards the new stable steady 
state. The associated time scale is approximately given 
by the transit time Tt r of the domain wall which connects 
the charged and the neutral regions. The general discus- 
sion of front propagation into unstable states [|l8| goes 
beyond the purpose of this paper. Here, we give only a 
rough estimate for a steep charge step 



T 2 

v \ V 



r 2 E 



(12) 



In particular, we neglected diffusion which acts to slow 
down the domain wall velocity and which acts to smear 
out the domain wall. Equation ( [12] ) can be obtained by a 
projection onto the translation mode of the domain wall 
and has the simple interpretation that the front travels 
with the drift velocity of the carriers. 
It should be noted that the slowness of the charge redis- 
tribution indicates a strong dependence on weak pertur- 
bations of the homogeneous insulator bulk. While weak 
forces are not expected to hinder the growth of the fast 



unstable injection mode, the charge redistribution can be 
considerably influenced by traps, grain boundaries etc. 
Therefore an experimental observation of the slow dy- 
namics and the final state discussed below requires a suf- 
ficiently clean material. Macroscopic insulating bodies 
used in high-voltage devices where the injection insta- 
bility occurs, are usually not very clean. Modeling the 
dynamics of the field limiting space charge should thus 
include bulk inhomogeneities. 

In Fig. we show a numerical simulation of charge in- 
jection in the perfect insulator for d = 1 and a — 3. After 
an increase of the voltage beyond instability threshold, 
a charge domain forms. On a long time scale the charge 
cloud smears out and relaxes finally towards a 1/r-like 
distribution, whereas the electric field becomes spatially 
uniform. In order to discuss this final steady state in the 
framework of our model, we first consider the case D = 0. 
From V • j = one finds a bulk solution 



E{r) = AdA 1 -^/^ 



p(r) = eA d 



l + da 
1 



-{a+d) / (l+a) 



a 



(13) 
(14) 



with Ai = V/(r 2 - n) and A 2 = V/(r^ /{1+a) - 
1+ )/(l + a^ 1 ) for the cylindrical and the spheri- 
cal case, respectively. Note the similarity of the electric 
field distributions for d = 1, a arbitrary, and for a — > oo, 
d arbitrary. In these cases, E(r) relaxes eventually to 
a constant value (see also Fig. Ha). Clearly, the bulk 
solution ( pl| ) does not satisfy the boundary conditions. 
For a small diffusion constant, the solution is expected 
to be changed considerably only in a boundary layer near 
the contacts. We find that the solution deep in the bulk 
far away from the contacts is only weakly disturbed by 
a small diffusion constant. For the cylindrical geometry 
(d = 1) and for V < r 2 E , the bulk solution reads in 
leading order of D 



E(r) 



r 2 avr V 



(15) 



For d = a = 1, Eq. (|l5| ) is the exact bulk solution. 
For a = 1 we do not find a linear instability of the ideal 
insulating state neither numerically nor with the analyti- 
cal approximation (|Io[). However, we conjecture bistabil- 
ity of neutral and charged state for V > (D/p,) \n(r 2 /ri) 
and a loss of global stability at a certain field. Below, 
this conjecture will be confirmed with the help of a sim- 
ulation. A detailed investigation of this case, however, 
will be published elsewhere. 

Assuming r 2 ^> n, Eqs. (|l^) and ( |l4| ) yield the current 



I d = b d evE 



V 

r 2 E 



a+l 



(16) 



where bi = 2nL z and b 2 = 4nr 2 a 1+a (l + 2a)/ (1 + a) 2+a 
for the cylindrical and the spherical case, respectively. 
Just above instability, the current is finite, though it is 
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small (of 0{{n/r 2 ) a+1 )). We mention that for d = 2, 
a = 1 and without diffusion, the l/\/r behavior of the 
electric field (13) is discussed by Lampert and collabo- 
rators H[[o|]. Neglecting diffusion, they obtain a bound- 
ary layer due to an E = boundary condition at the 
contact. One recovers from Eq. (^6|) in this case the 
current-voltage characteristic of the perfect insulator in 
a spherical conductor, / = (3/2)ire/j,V 2 /r2- A detailed 
discussion of the steady state, e.g. in the presence of 
intrinsic carriers, can be found also in Refs. P.tL0[. 



V. THE AC DRIVEN INSULATOR 

The localized field limiting space charge which forms 
at instability can be observed in ac experiments ||f|. 
Since A depends exponentially on a, the injection mode 
grows infinitely fast at E c = Eq in the limit of infinitely 
sharp mobility edge (a — * oo). Consequently, the electric 
field saturates immediately at E = E c due to sceening. 
On the other hand, the characteristic time of the charge 
redistribution, Eq. ([l^), diverges, provided E c ri > V. 
For E c r2 < V the whole bulk is charged up. 
In the following, we consider a periodic voltage which 
vanishes for t < and which, for positive t, is given 
by V(t) — Vsm(u>t). The frequency u> = 2ir/T obeys 
Lur tr ^> 1 ^> u)/X, where A is the steady-state stability 
eigenvalue for an electric field equal to the amplitude of 
the electric field oscillation. Since A -1 -C 1ms and the 
transit time r tr is of the order of hours [|J, reasonable 
frequencies are in the range of 1CP 1 — 10 3 Hz. 
Typical solutions are shown in Fig. |] for the cylindrical 
geometry and for various values of a. The thin solid line 
represents the reference electric field E\ /Eq at r± of the 
ideal insulating state with a purely capacitive response. 
The other curves represent numerical simulations of the 
ac response in the presence of injection. For fields be- 
low instability threshold, E\ < E c , the sample remains 
locally neutral. An increase of the field beyond thresh- 
old leads to the injection of charge in such a way that 
the local electric field at the inner contact is saturated 
slightly below the critical field. Negative and positive 
charge is periodically injected for a — 3 (dotted curve). 
Clearly, due to the electron-hole symmetry the solutions 
are symmetric with respect to inversion of the sign of 
the amplitude. The critical field is about 2Eq : where the 
field drops fast and saturates below Eq . This discontinu- 
ous transition from the neutral state to the charged state 
indicates also bistability. In the limit of a sharp mobil- 
ity edge (solid curve; a — 51), the neutral state decays 
immediately at E\ — Eq and the electric field oscillates 
between ±Eq, as expected. 

On the other hand, for a = 1 the insulating state p = is 
linearly stable even for large amplitudes V . For certain 
initial conditions or in the presence of additional current 
noise, however, we find also periodic solutions where only 
positive or only negative charge is injected. An exam- 



ple for positive charge injection is given by the dashed 
curve. Charge injection occurs at a large field amplitude 
E\ w %Eq. Once charge has been injected, a part of 
it remains in the sample and decreases the value of the 
electric field (dashed curve) compared to the chargeless 
case (thin solid curve). During the negative half-cycle, 
the electric stress is thus enhanced. Due to the presence 
of this charge, a further injection occurs at a much lower 
field in the second cycle of the oscillation. 
For initial conditions with reversed sign and V — > — V we 
find injection of charge with a different sign. We conclude 
that there are (at least) three different attractors indicat- 
ing bistability in the stationary case. Consequently, even 
in a system with electron-hole symmetry, rectification is 
possible due to dynamical symmetry breaking of charge 
injection. 



VI. CONCLUSION 

We have investigated charge injection in a macro- 
scopic and perfect insulator and for cylindrical and spher- 
ical geometries of the electrodes. The injecting metal- 
insulator contacts are modeled on a hydrodynamic level 
with boundary conditions for the charge density. We 
showed that, depending on the nonlinearity of the mo- 
bility, the ideal insulating state p = is unstable against 
a charge injection mode at a critical field E c , Former 
theories on charge injection assume either an intrinsic in- 
stability H or force injection directly by boundary con- 
ditions incompatible with a charge neutral state P,^0[. 
These theories cannot predict a finite critical field for 
the charge injection instability. For a macroscopic metal- 
insulator contact which is usually not well-defined, phe- 
nomenological boundary conditions to the charge density 
(which serves as an order-parameter field) is appropriate. 
Our theory predicts not only a critical injection field, E Cl 
but reproduces also the experimentally observed increase 
of E c with decreasing radius of the injecting contact Q. 
For a constant mobility (a = 1), the ideal insulating state 
is linearly stable. But numerical ac simulations which are 
adiabatic on the fast time scale show a decay of the neu- 
tral state to a charged state. From this we concluded 
bistability and a loss of global stability of the ideal insu- 
lating state. 

For a > 1, the time evolution from the insulating to 
the charged state occurs on two clearly separated time 
scales. On a short time scale, the injection mode grows 
at the instability and screens the electric field enhance- 
ment. A localized charge cloud forms near the contact. 
On a long time scale this charge redistributes over the 
whole sample. The localized field limiting space charge 
can be investigated with the help of an ac-bias which os- 
cillates with a characteristic time lying between the just 
mentioned time scales. 

Furthermore, we discussed the charged steady state for 
small diffusion constants. The current-voltage character- 
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istics of these solutions are determined by the bulk prop- 
erties and are thus equivalent to earlier results by Lam- 
pert and coworkers |^,|l0| in the perfect insulator limit. 
Future work should address the following problems. 
First, the transport model must be refined to include ad- 
ditional physical effects such as the influence of intrinsic 
carriers, boundary states, traps and impurity-band con- 
duction, surface potential decay, and (bi-)polarons. The 
effect of traps enters already by part via the diffusion 
constant. A small intrinsic carrier density is expected to 
increase the stability of the neutral state due to a finite 
dielectric relaxation mode. Furthermore, electron-hole 
symmetry is not very realistic and one expects rectifica- 
tion in ac experiments |19[] . 

Of interest is also the inclusion of heat transport and the 
influence of the temperature, and of mechanical stresses. 
Another important task is the determination of the pa- 
rameters appearing in the boundary conditions for the 
charge density from microscopic models. This is reason- 
able, however, only for physically well-defined contacts 
which is usually not the case for typical macroscopic 
metal-insulator contacts. 

The following interesting problem concerns the case of 
a constant mobility without linear instability. Injec- 
tion should then be associated with a nucleation of the 
charged phase at the contact [^o| , probably via a critical 
droplet with a the shape similar to the injection mode 
(|Tf|). Such an injection mechanism is also possible for 
bistability at a > 1 in the region where the ideal insulat- 
ing state is 'supersaturated'. 

Finally, in order to quantitatively compare theoretical 
with experimental results, one should investigate geome- 
tries different from cylindrical and spherical symmetry 
as, e.g., a tip-plate arrangement |f| f2f| , [22f . In contrast to 
the simple finite difference methods used in the simula- 
tions of this work, finite element methods are more ap- 
propriate to simulate charge injection in such real two- 
or three-dimensional geometries. 
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APPENDIX 

We calculate the relation between the critical value A c 
and the boundary-condition parameter k by solving Eq. 
(0) with the boundary condition (||) for Sp at A = 0. We 
can assume Sp > 0. Integration of Eq. ([?]) leads to 

A c 5 _ d s = M 

with a constant . The general solution of this equation 
is 



Sp=[A 1 -A 2 



r dr 

ex P(' 



ad — 1 



) exp(; 



f 



Applying the boundary conditions (g) to this function 
yields two linear equations for the constants A\ and 
A-}. The existence of a non-trivial solution requires 
the vanishing of the determinant associated with these 
equations. This condition can be written in the form 
an 2 + ok + c = with constants a, b, c depending on A. 
The solution for negative k defines the stability bound- 
ary plotted in Fig. Qa. In a similar way one calculates 
the dependence of E c on n which is shown in Fig. [l]b. 
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FIG. 1. a) Largest eigenvalue A of the stability problem as 
a function of the control parameter A (kxi = —0.5 d = 1, 
a — 3). b) Eigenfunctions of the stability problem. The solid 
curve represents the marginal charge-injection mode (A = 0) 
at A = A c . Modes are more localized at the inner contact as 
A increases. 

FIG. 2. a): Critical values A c as a function of — «ri, 
for cylindrical (solid) and spherical (dashed) geometries 
(fi/r2 = 0.01). Different curves with decreasing stability 
threshold belong to a — 1,3, and 15. b) Critical value of 
the electric field at the inner electrode as a function of the 
size of the injecting electrode (a = 3; solid: d = 1, dashed: 
d = 2). 

FIG. 3. Evolution of a) the electric field distribution and 
b) the charge density distribution beyond stability thresh- 
old. The localized injection mode (dotted curve) grows up to 
a certain amplitude. The domain wall moves into the bulk 
(dashed-dotted curves) until the final steady state with a uni- 
form field (solid curve) is reached (d= 1, a = 3). 

FIG. 4. Time dependence of the ac field at the inner elec- 
trode for a = 1 (dashed), a — 3 (dotted) and a = 51 (solid); 
vri/D — 0.1. Solid thin curve: field without injection. 
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